Domain Walls Motion and Resistivity in a Fully-Frustrated Josephson Array 
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It is identified numerically that the resistivity of a fully-frustrated Josephson-j unction array is 
due to motion of domain walls in vortex lattice rather than to motion of single vortices. 



A square two-dimensional periodic Josephson-junction 
array (JJA) in a uniform transverse magnetic field, with 
half a flux quantum per plaquette, is a realization of a 
fully frustrated (FF) XY model M . Its ground state is a 
checkerboard pattern of plaquettes with currents flowing 
clockwise and anti-clockwise [pl.The chirality, defined as 
the sum of currents around plaquette in a clockwise direc- 
tion, has antiferromagnetic order. Apart from Kosterlitz- 
Thouless (KT) S transition as in the unfrustrated case 
the FF XY model has an Ising-like transition of chiral or- 
dering ||], Numerical @ - ^ and analytical || work 
indicates that critical temperatures of KT and Ising tran- 
sitions are very close or equal. Mon and Teitel stud- 
ied current- voltage curves of the FF JJA using Langevin- 
dynamics simulation. They found that the temperature 
dependence of the exponent a of the current voltage char- 
acteristics (V ~ I a+1 ) does not follow the form expected 
from the results of equilibrium simulations of helicity 
modulus of FFXY model ||. In particular there were 
no sign of discontinuity in a at the critical temperature, 
while the discontinuity was observed in the helicity mod- 
ulus Q . Mon and Teitel proposed that this discrep- 
ancy is because apart from vortex-antivortex pairs (KT 
excitations) the system has excitations in the form of do- 
main walls in vortex lattice (Ising excitations), motion of 
which also contributes to resistivity. Disagreement be- 
tween the exponent of the current voltage characteristics 
of FF JJA and helicity modulus of FFXY model has been 
also observed experimentally, when in some cases [ pd| no 
sign of discontinuity in a was seen , and in others the 
sharp jump in a was observed Jl2[ , but it was substan- 
tially below T c , expected from simulations. In contrast 
for the unfrustrated case |Tg| ] the sharp jump in a was 
observed at the temperature equal to T c , found in simu- 
lations. Some experiments flij on superconducting wire 
networks in a magnetic field have reported exponential, 
instead of power-law, I—V curves. The exponential form 
is expected theoretically [[lO) if resistivity is due to motion 
of thermally activated domain-walls in vortex lattice. 

In this article it is shown, by numerical simulations, 
that domain-wall motion dominates the resistivity at 
least at low temperatures and small currents. 

In simulations the modification of Langevin dynamics 



method of Refs. jnj and jlB) proposed by Falo, Bishop, 
and Lomdahl [jl7| was used. All junctions are assumed to 
have the same critical currents Iq and to be shunted by 
resistances R and all superconducting nodes to have a ca- 
pacitance to the ground C. The resulting set of dynamic 
equations (which follows from Josephson equations and 
the charge conservation law) is: 
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Here 9 n is the phase, P n is the electrostatic potential 
and c/ n is the charge of the superconducting node on 
site n, A mn = I s f™ Adl is the integral of the vector 
potential from node n to m (sum of A mn around the 
plaquette is equal to the number of magnetic flux quan- 
tums through the plaquette) . Summation is over nearest 



neighbors only. I^ n is a thermal noise current with: 
< >= (2T/R)S(t - t), 



(4) 



with T the temperature. 

Instead of the uniform injection of Ref. injection 
from the superconducting bars at the edges of the array 
(see Fig. yj) was used in this work in order to imitate 
the experiment Jl2| , |15(] . Free boundary conditions in 
the transverse direction where used for the same reason. 
The bars were connected to each of the nodes at the 
edges of the array by Josephson junctions. The exter- 
nal current was injected into the left bar and extracted 
from the right bar. It was observed JTsJ that for the FF 
JJA injection from bars leads to unreasonably low zero- 
temperature critical current in comparison to uniform 
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injection. The reason for this is that FF plaquettes adja- 
cent to bars have only three Josephson-junctions leading 
to bigger phase differences in the ground state and there- 
fore lower critical currents. In this work plaquettes adja- 
cent to bars were assumed unfrustrated (experimentally 
this would correspond to this plaquettes having an area 
much smaller than plaquettes in the bulk) and the zero- 
temperature critical current was the same as in the case 
of the uniform injection [I c = 0.35), though lower than 
the theoretical value for an infinite array (7 C = y2 — 1) 

PL Call. 





FIG. 1. Josephson-j unction array with superconducting 
bars at the edges. The external current is injected into the 
left bar and extracted from the right. Voltage is measured 
between the bars. 



Simulations were done for the intermediate damping 
case, i.e. for McCumbcr parameter (3=1. The equa- 
tions of motion were integrated with discrete time steps 
At = 0.05 in units of inverse Josephson plasma fre- 
quency ujj = y/2elQ/TiC . Decreasing At ten times did 
not change the results. 
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FIG. 2. (a)Random staircase is the phase difference (in 
units of 2-7r) across the 16 x 16 unfrustrated array versus time 
(in units of inverse plasma frequency) for temperature T = 0.3 
and current I/Iq = 0.5. It is clearly seen that the phase dif- 
ference grows by 2n jumps, corresponding to vortex crossings 
shown in Fig. 4. (b)The same for the 17 x 17 fully-frustrated 
array for temperature T = 0.3 and current I/Io = 0.07. One 
can see the giant (~ 20 x 2n) jumps corresponding to domain 
wall propagation shown in Fig. 5. 



The phase difference cfi across the array was recorded 
as a function of time (see Fig. ||). The average voltage 
is: 



V 



^begin 



(5) 



where <t>begin arid (fiend are the phase differences across the 
array at the beginning and the end of the run correspond- 
ingly, and t run is the run time (in units of l/u>j). To 
estimate the errors each run was divided in four intervals 
and average voltages for this intervals were calculated. 
The estimated errors are -^f- ^0.1. 



In the following we shall adopt the following units: the 
current is measured in units of the critical current of the 
unfrustrated array Iq, the voltage in units of RIq, and the 
temperature in units of the Josephson energy hl{j/2e. 

Looking at Figure || (a) and (b), where phase difference 
across unfrustrated and frustrated arrays is plotted ver- 
sus time, one can notice a remarkable difference. While 
the phase difference across the unfrustrated array grows 
by 27r jumps, the one across the fully-frustrated array 
grows by giant ~ 20 x 27r jumps. In order to quan- 
tify results more we classify jumps according to their 
size. A jump of size 27rn is put into bin k such that 
2 fc_1 < n < 2 k . A fraction of the total phase slip, due 
to jumps in each bin is shown in Fig. We see that 
for unfrustrated array jumps in 1st and 2nd bins (that is 
n = 1, 2, 3) dominate, while for the fully-frustrated array 
jumps in 5th bin (n — 16, . . . , 31) give biggest contribu- 
tion to the total phase slip. 
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FIG. 3. A fraction, p, of the total phase slip in Fig. 2, 
due to jumps of size n2ir. Unfrustrated array - circles, 
fully-frustrated array - squares. 



As the Langevin dynamics simulation is performed us- 
ing pseudo-random numbers one can always go back and 
investigate what happens in the vicinity of the jump in 
detail. In Fig. 4 the 3 snapshots of the unfrustrated ar- 
ray in the region of one of the jumps are shown. One 
can clearly see a vortex crossing the array. Therefore 
the resistivity in the unfrustrated array is due to vortex 
crossing ^0|. In Fig. 5 the 3 snapshots of the FF array in 
the region of the giant jump are shown. It is hard to see 
what is going on in the FF array in the phase representa- 
tion, therefore the chirality (sum of currents around pla- 
quette in a clockwise direction) representation was used. 
However even then domain walls are not readily seen. To 
make them clearly visible one ground state is represented 
by vertical and other by horizontal arrows. One can see 
that before the giant jump the array is in one of chiral 
ground states, during the jump the domain of opposite 
chirality grows and at the end of the jump it fills all of 
the array. Resistivity of the FF array is therefore due to 
domain wall motion. 
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The method described here works, however, only for 
low temperatures (T < 0.35)and low currents (/ < 0.1). 
When the temperature or current is increased one has 
difficulties in identifying the separate jumps, while the 
phase difference - time curve is smooth. 



FIG. 4. Vortex crossing in the unfrustrated array of size 
16 x 16. Phases of superconducting grains are indicated. 
Three consecutive snapshots (a),(b),(c) are in the region of 
the phase jump at time t ~ 100 (see Fig. 2(a)). 
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FIG. 6. I-V curves (i = I/I vs v = V/RIo) of the 
fully-frustrated Josephson-junction array for 5 different tem- 
peratures (given in units of Josephson energy, hIo/2e). 
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FIG. 5. Domain wall propagation in the fully-frustrated ar- 
ray of size 17 x 17. Chirality is indicated. Two ground states 
are represented by vertical and horizontal arrows. Three con- 
secutive snapshots (a),(b),(c) are in the region of the giant 
phase jump at time t ~ 1.5 x 10 4 (see Fig. 2(b)). 



FIG. 7. Exponent a of the current-voltage characteristics 
(V ~ I a+1 ) versus inverse temperature (1/T), obtained from 
the data of Fig.6. Dashed line is obtained from helicity mod- 
ulus (Ref. [4]) using Eq. 6. 

The I — V curves of the FF array for different tem- 
peratures are shown in Fig. 6. And the exponent a 
(V ~ -^ a+1 ), obtained from them, is shown in Fig. 7. 
There is a good agreement between the data of Fig.6 and 
those of Mon and Teitel (Fig. 1(b) of Ref. goj). There is 
less agreement between the exponent a obtained in this 
work (Fig. 6) and that of Mon and Teitel (Fig. 2 of Ref. 
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pot). This is due to the fact that much lower currents 
have been studied in the present work than in Ref . [[l0| , 
where a was obtained by fitting the I — V curves in the 
region of current just below the critical one. In this work 
a was obtained by fitting in the region of current 0.05-0.1 
(in units of the critical current of the unfrustrated array). 
For higher current single-junction effects are important, 
for lower ones - finite size effects ^lj cause the change 
of the I — V curve from nonlinear to linear. It is inter- 
esting, that though the resistivity is due to domain walls 
motion the I — V curves still can be fitted by the power 
law. In Fig. 7 the theoretically expected [£2| exponent 
of the I — V relation, calculated as 



a = 7rr/T, 



(6) 



using helicity modulus T obtained in Ref. ^ is shown 
by dashed line. Though agreement between a predicted 
using helicity modulus and obtained from I — V curves of 
the present work is better than in Ref. [^0) , still there is 
disagreement. The reason for it is probably domain wall 
contribution to resistivity, as suggested by Mon and Tei- 
tcl jnj and demonstrated explicitly in the present work. 

In conclusion in this work an approach is suggested 
when one records phase difference across the Josephson 
array versus time (time derivative of the phase difference 
is the voltage). Resulting pictures are different in the 
cases of unfrustrated and fully frustrated arrays. In the 
first case phase difference grows by 2n jumps, while in the 
last by giant (many 2n) jumps. Snapshots of the array 
in the region of a giant jump show that during it there 
is a transition from one of two chiral ground state to an- 
other, therefore identifying the mechanism of resistivity 
as domain walls motion. 
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